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Abstract 

Transmit signal and bandwidth optimization is considered in multiple-antenna relay channels. Assuming all 
terminals have channel state information, the cut-set capacity upper bound and decode-and-forward rate under 
■ full-duplex relaying are evaluated by formulating them as convex optimization problems. For half-duplex relays, 

I bandwidth allocation and transmit signals are optimized jointly. Moreover, achievable rates based on the compress- 

^ ■ and-forward transmission strategy are presented using rate-distortion and Wyner-Ziv compression schemes. It is 

I observed that when the relay is close to the source, decode-and-forward is almost optimal, whereas compress-and- 

forward achieves good performance when the relay is close to the destination. 

00 _ 

Index Terms 

Bandwidth allocation, convex optimization, relay channels, multiple-antenna, transmit covariance matrices. 

^ ■ I. Introduction 

In wireless communications, the transmission rate is fundamentally limited by the channel propagation loss over 
the range of transmission. Relaying has been proposed as a transmission strategy that can improve the performance 
^ of wireless systems. In a relay channel, in addition to the source and the destination, there is also a relay terminal. 
^ The relay does not have its own data to send or receive; its intention is to facilitate the transmission between 
Q\ the source and destination. On one hand, the relay brings additional power to the network, as typically, the relay 
CN is under its own power source. On the other hand, a relay can also help to shorten the transmission range by 
enabling communication over two hops, and cooperate with the source to perform joint encoding of the transmit 
signals. However, if the relay transmission scheme is not designed properly, the relay may also create undesirable 
interference to the terminal at the destination. In this paper, we investigate the optimization of transmit signals 
in relay channels. In particular, we consider a multiple-input multiple-output (MIMO) relay channel, where the 
• 'H source, relay, and destination terminal all have multiple antennas. In wireless communications, using multiple 
rS transmit and receive antennas has been shown to provide substantial improvement in channel capacity Ol, We 
d also evaluate the MIMO relay channel rates under different network geometry, and investigate the effectiveness of 
the corresponding relaying schemes. 

The three-node relay channel model is proposed in |3|. In |4|, a capacity upper bound and achievable coding 
strategies are presented for the relay channel, but the relay channel capacity remains an open problem. For Gaussian 
single-antenna relay channels, capacity bounds and power allocation are studied in [5|. Capacity bounds on half- 
duplex relaying are presented in |61, |[7|. Bandwidth and power allocation are considered in [81, [9J for fading 
orthogonal relay channels, and in |[TOl for the amplify-and-forward scheme in Gaussian relay networks. Transmitter 
cooperation versus receiver cooperation in relay channels is compared in lITTI . Relay channel coding strategies, with 
extensions to relay channels with multiple terminals, are given in lfT2l . For relay channels with multiple-antenna 
terminals, bounds to the cut-set capacity upper bound and decode-and-forward rate are considered in fl3l . 114 1. In 
iflS I. [16 1, the diversity-multiplexing tradeoff is characterized for full-duplex and half-duplex MIMO relay channels. 

In this paper, we consider a multiple-antenna relay channel where the terminals have knowledge of the channel 
state information (CSI). We consider optimization of the source and relay transmit signals to evaluate the MIMO cut- 
set capacity upper bound and the decode-and-forward achievable rate, by formulating them as convex optimization 
problems. In the case of half-duplex relaying, the bandwidth allocation and the multiple-antenna transmit signals 
are optimized jointly. We also present achievable rates in MIMO relay channels using the compress-and-forward 
approach, under which the rate-distortion and Wyner-Ziv compression schemes are considered. 
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Fig. 1. Multiple- antenna relay channel. 



The remainder of this paper is organized as follows. Section JI] presents the multiple-antenna relay channel 
model, and the capacity upper and lower bounds. The cut-set bound and decode-and-forward rate optimization 
formulations are described in Section JII] under full-duplex assumptions, while Section |IV] considers half-duplex 
relaying where the relay cannot simultaneously transmit and receive in the same frequency band. The compress- 
and-forward transmission strategy is studied in Section |Vl and Section |Vl] concludes the paper. 

Notation: In this paper, H (]R,+, is the set of real numbers (nonnegative, positive real numbers), C is the 

complex field, 1 denotes the two-element set {0, 1}. Dimensions of vectors/matrices are indicated by superscripts. 

is the set of x positive semidefinite Hermitian matrices. X )^ Y {X '^Y) means the matrix X — Y i& 
positive (semi)definite. /at is the N x N identity matrix. and are the transpose and conjugate transpose, 
respectively, of a matrix A. The operators E[-], det, tr denote, respectively, expectation, determinant and trace. 
For random variables, x CM{n, Q), where x,fi£€^,Q £ M^, means that x is a circularly symmetric complex 
Gaussian random A^-vector about mean /u with covariance matrix Q. 

II. System Model 

A. Channel Model 

Consider a three-node wireless relay channel as illustrated in Fig. [T] The source node wishes to send a message 
to the destination; the relay node does not have its own message to send, but facilitates the transmission between the 
source and destination. Suppose the source has Mi transmit antennas, and the destination has A^i receive antennas. 
We assume the relay has M2 transmit antennas and receive antennas (for instance, in full-duplex operation the 
relay may have different sets of transmit and receive antennas). We consider a discrete-time flat-fading channel 
model, which is described by 

yi = Hiixi + H12X2 + zi (1) 

y2 = H21X1 + Z2 (2) 

where xi E C^^\ X2 € C^^^ are the respective transmit signals of the source and relay; yi G y2 G are 
the respective receive signals of the destination and relay; and zi ~ C7V(0, /atJ € Z2 CM{0,I]\[2) € C^^ 
are independent zero-mean circularly symmetric complex Gaussian (ZMCSCG) noise at the destination and relay, 
respectively. The complex baseband channel from the source to destination is Hn G c^ixA^i; from the source to 
relay is H21 G C^^^^^^; and from the relay to destination is H12 G C^i^*^^ 

We consider a block-fading channel model: the channels realize independently according to their distribution at 
the beginning of each fading block, and they remain unchanged within the duration of the fading block. In this 
paper, we assume the channel states can be estimated accurately and conveyed timely to all terminals: i.e., we 
assume channel state information (CSI) is available at all nodes. The source and the relay are under the respective 
transmit power constraints: E[3;^3;j] < Pi, i = 1,2, where the expectations are over repeated channel uses within 
every fading block. Power allocation across fading blocks is not permitted. The transmit signals have zero mean: 
E[xi] = 0, E[x2] = 0. 

It is convenient to write ([T]), ^ in block-matrix form 



y = Hx + z 



(3) 
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where 



H 



Hii Hi2 
H21 



G C 



iVxAf 



Xl 
X2 



G C 



N 



(4) 



with M = Ml + M2, and = A'^i + Moreover, we denote the joint covariance matrix of the transmit signals 
of the source and relay as 



Qu Q12 
Q21 Q22 



Q = E[xx^] G H 



M 



(5) 



where the conformally partitioned blocks (with respect to xi,X2) of Q have dimensions 

Qn = E[xixf ] G Q22 = E[X2X^] G Uf, Qu = E[xix^] = Qg G (6) 



B. Capacity Bounds and Achievable Rates 

The capacity of a relay channel is, in general, an open problem; however, there are known upper and lower 
bounds. The cut-set bound described in lITTl provides an upper bound to the relay channel capacity. Intuitively, 
the cut-set bound states that, over any possible joint source-relay transmit signals, the relay channel capacity cannot 
exceed the smaller of i) the maximum rate at which information can flow out of the source, and ii) the maximum 
rate at which information can flow into the destination. 

On the other hand, capacity lower bounds of the relay channel, as achieved by two coding strategies, are given in 
0. In the decode-and-forward strategy |T, Thm. 1], transmission is done in blocks: the relay first fully decodes the 
message from the source in one block, then in the ensuing block, the relay and the source cooperatively transmit 
the message to the destination. In the compress-and-forward strategy \4. Thm. 6], the relay does not decode the 
source's message but sends a compressed version of its observed signal to the destination. The destination then 
combines the compressed signal with its own receive signal to decode the source's message. A detailed discussion 
on the relay channel coding strategies can be found in [12|. 

In the following sections, we evaluate these capacity bounds and achievable rates for full- and half-duplex 
Gaussian MIMO relay channels. Under half-duplex relaying in Section |IVl we impose the constraint that the relay 
cannot simultaneously transmit and receive in the same frequency band. For performance comparisons, we also 
consider: direct transmission when the relay is not available (i.e., the capacity of an Mi x Ni MIMO channel 
under transmit power constraint Pi); orthogonal two-hop relaying (Section [iV-CI l; and the scenario when the relay 
is co-located with the source or destination (Section IIV-DI ). 



III. MIMO Relay Channel Capacity Bounds 

In this section, we present the optimization frameworks for evaluating the MIMO relay channel cut-set capacity 
upper bound and the decode-and-forward achievable rate. We first adopt the full-duplex assumption where the relay 
can transmit and receive in the same frequency band at the same time. In practice, full-duplex transmission is difficult 
to realize. Nevertheless, the full-duplex model provides insight into the design of effective coding strategies for relay 
channels, and serves as a performance upper bound for half-duplex systems. Half-duplex relaying is considered in 
Section |IV] where the relay cannot simultaneously transmit and receive in the same band. 



A. Cut-Set Capacity Upper Bound 

The cut-set capacity upper bound 14] for the relay channel is given as an optimization in terms of the channel 
mutual information as follows: 

Rcs= max mm{l{xi;yi,y2\x2), Iixi,X2;yi)} (7) 

P{Xl,X2) 

= max min{/(xi; 2/1x2), (8) 

p{x) 
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Gaussian signals are optimal in the cut-set bound and decode-and -forward rate |[T2l Proposition 2]. We denote the 
transmit signals by: x ~ CA/^(0, Q), where Q G H^^ is the covariance matrix of x. Then the mutual information 
expressions in ([S]) evaluate |2] to 



Rcs= max min{logdet(/iv + ^i<5i|2^^f ), logdet(/jVi )} 

Q : tr Qii<Pi, i=l,2 

where Hi and are, respectively, the first block column and block row of H 



Hi 



Hn 
H21 



G C 



NxMi 



Hi ^ [Hn Hu] e C 



NixM 



(9) 



(10) 



and the conditional covariance matrix Qi\2 — E[xix(^|x2] is given by the Schur complement of Q22 in Q 

Qi\2 = Qii - Qi2Q22Q'ii (11) 

where we assume Q22 >~ 0. The zero-mean Gaussian signal x is fully characterized by its covariance; therefore, in 
Q, the sole optimization variable is the joint covariance matrix Q. 

The cut-set bound maximization in Q can be formulated as the following optimization problem: 



maximize i?cs 

over i?csGlR+, QeHf, Qi|2GHf^ 
subject to i^cs < logdet(/Ar + HiQi\2H^) 

i?cs <logdet(/iv, +^iQ^f) 
tr(Cf gCi) < Pi 
tv{C^QC2) < P2 

Q - C1Q112CI h 



where Ci , C2 are constant matrices defined as 



Ci 







G 1 



MxMi 



Co 





hu 



G 1 



MxMn 



(12) 
(13) 
(14) 
(15) 
(16) 
(17) 
(18) 

(19) 



In the optimization, ([141 ). ([TSl ) follow from the two terms inside the min expression in (|9ll; and ([T6l ). ([TT] ) represent 
the per-node transmit power constraints at the source and relay, respectively. The constraint ([TSl l results from 
relaxing the equality constraint in (fTTl) 

|2 = Qll — Ql2Q22^Q21 =^ <5l|2 ^ Qll — Ql2<522^Q21- 



Ql\1 



(20) 



By the semidefiniteness property of Schur complements |[T8l . we have the following identity on the right-hand side 
of (EOl): 



{Qll — Ql\2) — Ql2<522^Q21 ^ 



(Q11-Q1I2) 
Q21 



Q12 
Q22 



^ 



(21) 



where the right-hand side of (|2TI ) is equivalent to Q — CiQi\2Cf >: when written in the block-matrix form as 
defined in ([5]l. Finally, we show the relaxation in (|20l ) does not increase the optimal value in (fT2l ). Suppose given 
a set of fixed Qn, Q12, Q2i,Q22, we consider all X G H:Jl^^ such that X ^ Qii — Qi2Q22Q2i- Recalling that the 
determinant is matrix increasing |fT9l on the set of positive semidefinite matrices, we get 

log det(/^ + HiXH^) < log det(/^ + Hi{Qu - Qi2Q22Q2i)H^) (22) 

which only limits the feasible set in ([141) . 

The maximization in (fT2l)-([T8l) is a convex optimization problem; in particular, the log-determinant function is 
concave on positive definite matrices |[T9l . The solution of (fT2l ) can be efficiently computed using standard convex 
optimization numerical techniques, for instance, by the interior-point method |[T9l . ||2p^l. The above optimization 
can also be solved by the CVX 11211 . ll22l software package, which uses a successive approximation approach to 
model the log-determinant inequalities. All optimization formulations presented in this paper are convex problems, 
unless otherwise noted. 
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In addition to the per-node power constraints ( [T6b . ( fTTl ). if the source and relay are also under per-antenna 
power constraints, they can be readily incorporated in the convex optimization formulation. Let the antenna power 
constraints of the source and relay, respectively, be • • • and P2,i, ■ ■ ■ ,Pi,M2- The per-antenna power 

constraints are represented by 

qi,i < Pi,i, i = 1, . . . , Ml, qM,+j,Mi+3 < j = 1, . . . , M2 (23) 

where qi^i is the (i, i) entry of the co variance matrix Q, with i = 1, . . . , M. 



B. Decode-and-Forward Achievable Rate 

The decode-and-forward lH Thm. 1] relay channel achievable rate is given by 

i?DF = max mm.[l{xi]y2\x2), I{xi,X2;yi)] (24) 
= max min{/(xi;y2|a:;2), (25) 

p{x) 

max min{logdet(/jv, +-ff2iQi|2-ff|f), logdet(/7Vi +^iQ-H'f )} (26) 

Q : trQii<Pi, i=l,2 

where Qi\2 is as given in ([iTI ). The decode-and-forward rate can be formulated as the solution to the following 
convex optimization problem: 



maximize Rdf (27) 

over i?DF e Q e Hf , Qi|2 G ^ (28) 

subject to R-DF < logdet(Iiv, + iJ2iQi|2^^K) (29) 

i?DF < log det(/^, + HiQH^) (30) 

tr(Cf QCi) < Pi (31) 

tr(CjQC72) < P2 (32) 

Q - CiQi|2Cf ^ 0. (33) 

The derivations for the optimization problem formulation parallel those presented in the previous section. Note 
that the decode-and-forward optimization problem (|27])-([33]) is similar to the cut-set bound optimization problem 
(fT2l)-([T8]). The only difference is the aggregate channel Hi in (fT4l ) versus the source -relay channel H21 in ( |29l ). 
Consequently, we expect the decode-and-forward scheme is nearly capacity-achieving when the direct channel Hu 
is weak relative to the source-relay channel H21. 



C. Numerical Results 

In the numerical examples in this paper, we assume a network geometry as depicted in Fig. |2l In the two- 
dimensional network, the source is located at coordinates (0,0), the destination is at (1,0), and the relay is at 
{dx,dy). We will use a distance-based path-loss power attenuation exponent rj = A, combined with independent and 
identically distributed (i.i.d.) Rayleigh fading for each channel matrix entry. Shadow fading can also be included, 
but it is omitted here to allow a simple geometric interpretation of the relay network topology. The relay channel 
matrices are given as 

Hn = H(J\ H2, = (dl + dlr^/'H(J\ H12 = ((1 - d^f + 4)-''/'/7(f) (34) 

where H^^^ G ^iVixMi^ ^^2) ^ ^n,xM,^ ^ ^n,xM,^ ^^^^ ^^^^ ^^^^ of j^W ^ ^^2)^ ^(3) _ ^^^^^ 

For the numerical experiments presented in this paper, 50 random instances of the channel realizations Hw \ Hw\ 
are generated. Then under each channel realization, the corresponding optimization problems are solved to 
evaluate the relay channel capacity bounds and achievable rates. For consistent comparison, the same set of channel 
realizations is used to compute the performance of the different coding schemes under consideration. The convex 
optimization problems are solved using the barrier interior-point algorithm described in [|19i Section 11.3]. 

The empirical cumulative distribution functions (CDFs) of the cut-set (CS) bound and decode-and-forward (DF) 
rate of a MIMO relay channel, where the relay is located at {dx^dy) = (1/3, 1/2), are shown in Fig. [3](soUd lines). 
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Fig. 2. Geometry of the source, relay, and destination nodes. 




Fig. 3. Full-duplex CDFs of MIMO relay channel cut-set (CS) bound and decode-and-forward (DF) rate. The relay is located at (1/3, 1/2). 
All terminals have four antennas: A'li — Ni — M2 = N2 ~ 4, and unit power constraints: Pi = P2 = dB. The dotted lines represent the 
rates under per-antenna power constraints. The plots (a) and (b) correspond to the upper and lower bounds, respectively, from 1131 Thms. 3.1 
and 3.2]. 



All terminals have four antennas: Mi = Ni = M2 = N2 = 4:, and unit power constraints: Pi = P2 = dB. 
The dotted Unes represent the rates under per-antenna power constraints, where all antennas have equal constraints: 
pi^i = ■ ■ ■ = pi A/i = Pi/Mi, and p2,i = • • • = P2,M2 = P2/M2. Also shown in the plot is the MIMO capacity 
of the direct channel Hu when the relay node is not available (No Relay). The decode-and-forward achievable 
rate considerably outperforms the direct channel capacity and is quite close to the cut-set capacity upper bound. 
For all coding schemes, imposing per-antenna power constraints only slightly reduces the rates as compared to 
per-node power constraints. For comparison, the upper and lower bounds from |13, Thms. 3.1 and 3.2] are plotted 
and labeled (a) and (b), respectively. It is observed that the capacity upper and lower bounds can be tightened when 
the transmit signals of the source and relay are optimized. The upper bound (a) is computed by searching over 
combinations of p = {0, 0.05, . . . , 0.9, 0.95} and a = {10^^, 10^*^'^, . . . , 10°'^, 10^}, where the parameters p, a are 
as defined in fT3ll . For each choice of p,a, a convex optimization problem is solved. 

IV. Half-Duplex Relaying 

In Section |llll we assumed the relay was able to transmit and receive simultaneously in the same band. Such 
full-duplex radios can be difficult to implement in practice. In this section, we consider a half-duplex relay, where 
the relay receives in one band and transmits over a different band. In particular, we assume the channel has unit 
bandwidth, and it is partitioned into sub-channel Band 1 with bandwidth wi, and another orthogonal sub-channel 
Band 2 with bandwidth W2, with wi + W2 < 1. The relay can only receive in Band 1 and it can only transmit in 
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Band 2. Hence the channel is described by 

_ TT (1) , ^(1) „{2) _ rr (2) , rj (2) , (2) 

Ui — J^uXi + , — Hiix^ + H12X2 + (33) 

_ TT Jl) , (2) (2) _ , . 

where the superscripts designate the corresponding bands. The noise powers in the sub-channels are given by 

E[zi'\zi'^f] = wiIn., E[zf )(zf ))^] = W2lN„ i = 1,2. (37) 

Let Q^ii , Q^'^'^ be the transmit signal covariance matrices in the two bands 

We assume the bandwidth allocation and the transmit signal covariances in each band can be optimized with 
respect to the channel realizations. In the following sections, we consider the cut-set capacity upper bound and the 
achievable rates under the half-duplex relaying constraint. 



A. Half-Duplex Cut-Set Bound 

Let the mutual information across the different cut sets be designated as labeled in Fig. |4l The cut set around 



the source is shown in Fig. |4(a)| Let Ri,R2, respectively, denote the egress information rate out of the source in 
Band 1 and Band 2. On the other hand, for the cut set around the destination shown in Fig. |4(b)[ let Rd,Rc be 
the ingress information rate into the destination in Band 1 and Band 2, respectively. Optimizing over the transmit 
signals and the bandwidth allocation, the half-duplex cut-set bound is characterized as follows: 

maximize i?hcs (39) 

over Ri,cs,Ri,R2,Rd,Rc,wi,W2 € 1R+, qS^^ € H^S Q^^) ^ jgM (4^^ 

subject to i?hcs < min(i?i + R2, Rd + -Rc) (41) 

Ri < wi log det {In + ir.HiQiiH^) (42) 

R2 < W2 log det(/,v, + ^^HnCfQ^^^CiHii) (43) 

Rd < wi log det {In, + ^r^ii^ff ^n) (44) 

Rc < W2 log det(/jv, + ^^HiQ^^^Hi") (45) 

tr gff + tr(Cf < (4g) 

tr{C^Q^^^C2) < P2 (47) 
wi + W2<l (48) 

where Ci,C2 are as defined in ([T9l ). By continuity we define: ?i;logdet(/ + X/w)\w=q = 0, for all X ^ 0. The 
right-hand side of each constraint in (l42l ). (1431 ). (l44b . (1451 ) is a concave function, being the perspective of the log- 
determinant function. (Given a function /(x), the perspective of / is defined as the function g{x,t) = tf{x/t), 
t G 11++, and the perspective operation preserves convexity |[T9l .) 



B. Half-Duplex Decode-and-Forward Rate 

Fig. |5] depicts the operation of decode-and-forward in the half-duplex mode. In Band 1, the source sends to the 
relay at rate R,-, of which R^ is decodable at the destination. The relay fully decodes the message from the source, 
and in Band 2 the source and relay cooperatively send to the destination additional information at rate R^. The 



8 




Source Destination Source Destination 

(a) Cut set around the source. (b) Cut set around the destination. 



Fig. 4. Half-duplex cut-set bound. 




Source Destination 



Fig. 5. Half-duplex decode-and-forward. 



half-duplex decode-and-forward optimization is given as follows: 

maximize i?hDF (49) 

over i?hDF, ^r, Rd, Rc, wi,W2€ R+, Q^u' e , Q^^^ e (50) 

subject to i?hDF < min(i?r, + Rc) (51) 

Rr < log det {In, + J-i/siQ^i/lf) (52) 

Rd < wi log det(/iv, + ^i^iiQff i^il) (53) 

Rc < W2 log det(/^, + J-i7iQ(2)#f ) (54) 

tr Q^^l + tr(Cf < ^55^ 

tr(CjQ(2)c2) < P2 (56) 

u'i+u'2<l. (57) 



C. Two-Hop Relaying 

Two-hop relaying is a simple scheme that imposes relatively small coordination overhead between the source and 
relay. Its operation is portrayed in Fig. |6] In Band 1, the source transmits to the relay with signal covariance Qn''- 
The relay then decodes the message from the source, and re-encodes it to transmit to the destination in Band 2 



with covariance Q22 £ The following rate is achievable 

maximize i?2hop (58) 

over i?2hop, i?sr, R,d, wi,W2e R+, qSi^ G ^ , G Hf = (59) 

subject to i?2hop < min(-Rsr, ^rd) (60) 

R,r < wi log det(/jv, + ^H2iQ^Sh2i) (61) 

Rrd < W2 log det [In, + ^^i2Qg^ H^2) (62) 
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Fig. 6. Two-hop relaying. 



trQSl^<Pi (63) 

trQg<P2 (64) 

u;i + u;2 < 1 (65) 

where i?sr represents the rate from the source to the relay, and R^d represents the rate from the relay to the 
destination. 



D. Co-Location MIMO Capacity 

For comparison, we also consider the performance of the relay channel when the relay is close to the source 
or the destination. When the relay is co-located with the source, where it can cooperate perfectly with the source 
without overhead, the resulting channel is equivalent to an M x A^i MIMO channel. In the MIMO channel. Mi of 
the transmit antennas are under a sum power constraint of Pi, and M2 antennas are under sum power constraint 



P2. The MIMO capacity is given by 

maximize RmxNi (66) 

over RMy.N, € 1R+, Q € Hf (67) 

subject to i?MxAfi < logdet(/Ar^ +#iQFf ) (68) 

tr{CjQCi) < Pi (69) 

tr{C^QC2) < P2- (70) 



The above maximization is a convex optimization problem, and can be solved by the software package SDPT3 
II23I . which directly supports the log-determinant construct in the optimization objective function. 

On the other hand, when the relay is co-located with the destination, we assume they can cooperate perfectly 
without overhead. In this case, the resulting channel is equivalent to an Mi x N MIMO channel under a transmit 
power constraint of Pi. The MIMO channel capacity is 

RM,y.N= max logdet(/jv + //iQii-H'f) (71) 

trgii<Pi 

where the solution is given by waterfiUing power allocation ll24l along the eigenmodes of H^Hi. 



E. Numerical Results 

Fig. |7] shows the empirical CDF of the half-duplex cut-set (hCS) bounds, the half-duplex decode-and-forward 
(hDF) rates, and the two-hop relaying (2hop) rates for the MIMO relay channel with parameters as described in 
Section IIII-CI Again, the dotted lines represent the rates under per-antenna power constraints. For comparison. 
Fig. |7] also includes the full-duplex cut-set bound, full-duplex decode-and-forward rate, and the MIMO capacity of 
the direct channel Hu without the relay. It is observed that the half-duplex rates fall moderately as compared to 
the full-duplex rates. Moreover, the gap between the decode-and-forward rate and the cut-set capacity upper bound 
widens under the half-duplex mode. On the other hand, half-duplex decode-and-forward still provides a sizable 
capacity gain over direct transmission, while the two-hop relaying scheme achieves only marginally higher rates 
than when the relay is not available. Similar to the full-duplex case, imposing the per-antenna power constraints 
reduces the rates only slightly. 

Next, we investigate the relay channel capacity bound and achievable rate as a function of the relay position. 
In the following numerical experiments, we fix dy = 1/10, and vary dx from -1/2 to II/2; therefore, the relay 
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Fig. 7. Half-duplex CDFs of MIMO relay channel cut-set (hCS) bound, decode-and-forward (hDF) rate, and two-hop (2hop) relaying rate. 
The relay is located at (1/3, 1/2). All terminals have four antennas: Mi = Ni = Ah — N2 = and unit power constraints: Pi = P2 = dB. 
The dotted lines represent the rates under per-antenna power constraints. 

ranges from being closer to the source, to being closer to the destination. Again, all terminals in the network have 
four antennas: Mi = Ni = M2 = N2 = 4:, and unit power constraints: Pi = P2 = dB. The average rates for 
the different full- and half-duplex relaying schemes are plotted in Fig. [H they are computed over the 50 sets of 
random channel realizations with distance-based path-loss as given in (l34l ). The 8x4 MIMO capacity given by 
(|66l ), corresponding to the case where the relay is co-located with the source, is indicated by a circle, whereas 
the 4x8 MIMO capacity (TtT] ). under relay-destination co-location, is indicated by a square. For the half-duplex 
schemes, the bandwidth allocation in terms of the relay location is shown in Fig. |9] 

Under the full-duplex mode of operation, the decode-and-forward relaying scheme offers substantial capacity 
gain over transmission using only the direct channel. Over a wide range when the relay is close to the source, 
the decode-and-forward rate almost coincides with the cut-set capacity upper bound, and it is close to the 8x4 
MIMO capacity when the relay is at {dx,dy) = (0,1/10). The highest decode-and-forward rate is attained when 
the relay is located approximately midway between the source and destination. However, as the relay moves from 
the source and approaches the destination, the decode-and-forward rate begins to deteriorate. In fact, when the 
relay enters into proximity of the destination, the decode-and-forward rate underperforms direct transmission. This 
is because the decode-and-forward scheme requires the relay to fully decode the message from the source, and 
consequently the source-relay channel becomes a performance bottleneck. In practice, the source would typically 
enlist the relay's help only if it offers a capacity gain over direct transmission; hence the achievable rate may be 
taken as the maximum of the relaying rate and the direct rate. 

The half-duplex decode-and-forward scheme follows similar trends, but exhibits a wider gap from the half-duplex 
cut-set capacity upper bound, and its maximum capacity gain over the direct channel is less pronounced. In Fig. |9j 
it is observed that the system bandwidth is disproportionately allocated to Band 1 when the relay is in the vicinity 
of the destination, which corroborates with the system performance being limited by the source-relay link. The 
orthogonal two-hop relaying scheme, on the other hand, does not perform as well as decode-and-forward. It only 
offers moderate capacity gain over direct transmission, where the gain similarly is at its peak when the relay is 
about equidistant from the source and destination. 




11 



14 
12 



CO 



CP 
fcuO 

<1 



-8^ 









- 


/ Vir"Q \\ 


- 


8x4 




- 




/^hop \\\ 


A O 

4X5 


No relay 



















0.5 

Relay location: 



1.5 



Fig. 8. Full- and half-duplex cut-set (CS) bound, decode-and-forward (DF) rate, and two-hop (2hop) relaying rate with respect to the relay 
location. 



^0.8 



0.6 

Id 



PI 

CD 
bjO 
oj 
u 

< 0.2 



0.4 



-8 









\hDF // / 








2hop ^ 









0.5 

Relay location: dx 



1.5 



Fig. 9. Bandwidth allocation in the half-duplex relaying schemes with respect to the relay location. 



12 



V. COMPRESS-AND-FORWARD RELAYING 

From discussion in the previous section, it is clear that requiring the relay to decode the source's message can 
become a performance bottleneck when the source-relay channel is weak. In this section, we consider the compress- 
and-forward [4, Thm. 6] strategy where the relay does not attempt to decode the message from the source. Rather, 
the relay forwards a compressed version of its observation to the destination. The relay's observation is compressed 
in the sense that a finite number of bits is used to represent the analog signal (the scheme is sometimes also referred 
to as quantize-and-forward). Unlike the cut-set bound and decode-and-forward formulations, however, the transmit 
signal design and bandwidth allocation under compress-and-forward do not appear to be convex problems. In this 
section, we consider achievable compress-and-forward transmission schemes. We focus on full-duplex transmission; 
under fixed bandwidth allocation, the compression-and-forward operation readily extends to half-duplex relaying. 

A. Compress-and-Forward Transmission 

We first describe the general compress-and-forward strategy; specific compression schemes are considered in 
Sections |V-B[ IV-Cl The optimal joint design of the transmit signals and compression rate appears to be intractable; 
in the following we present suboptimal approaches to consider specific power allocation and compression schemes. 
We assume the source and the relay use Gaussian signals. Using the capacity-achieving strategy as in a multiple- 
access channel ifTTl . suppose the destination performs successive interference cancellation to allow simultaneous 
transmission from the source and relay. In particular, we consider the decode order in which the destination first 
decodes the relay's message, treating transmission from the source as noise. Then the relay's codeword is subtracted 
from the observed signal, and the message from the source is decoded. The source-destination transmission is thus 
interference-free from the relay's signals, and the source optimizes its own transmit signal covariance Qn according 
to 

Rii= max \ogdet(lN,+ HiiQiiH^A (72) 

where the solution is given by the waterfilling procedure. Let denote the covariance matrix that maximizes 
(1721 ). Next, the relay optimizes its transmit signal against the interference from the source's transmission 

Ri2= max logdet(/jv, +^i2Q22^i1) (73) 

where H12 is the effective channel from the relay to destination treating interference from the source as noise 

^12 = {In, + HuQuHfiy'/'Hi2. (74) 

Similarly, the solution in (1731 ) is given by waterfilling against the effective channel Hi2. 

In the compress-and-forward approach, the relay sends y2 € to the destination, which is a compressed 
version of the relay's receive signal y2, with compression rate R12 as given in (TTSl l. The compression schemes 
considered in this paper can be modeled by 

§2 = Ay2 + z (75) 

where A € C^^^^^ is a constant scaling matrix, and z ~ CJ\f{0,Z) £ is independent additive Gaussian 
compression noise, with Z £ Upon receiving y2 at the destination, the relay network is equivalent to an 

Ml X (A^i+A''2) MIMO channel, except that of its receive antennas are scaled by A and corrupted by compression 
noise z 

= AH2,\ + [az2 + z\ ■ ^'^^ 
The corresponding channel capacity is given by 

RcF = log det {In + ^iQn^f ) (77) 

where Hi is the effective source-to-destination-and-relay MIMO channel, incorporating the degradation introduced 
by the compression scheme as given in (1761) 

{Z + AA^)-'/'AH2i 
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The different compression schemes considered in this paper differ in their respective achieved values of A and Z. 
They are described in more detail in the next sections. 

B. Rate— Distortion Compression 

Let §2 G represent the compressed version of the signal y2 under rate-distortion theory |fT7|. Let z = y2 — y2 
be the compression error, where z has zero mean. The quality of the compression is characterized by the distortion 
covariance matrix: D £ H^^ = E[zz^]. In general, we wish D to be small to obtain a high compress-and-forward 
rate. In the following, we consider a simple approach to model this dependency by considering minimizing: tr D, 
which represents the sum of the quadratic distortion measures in the entries of 2/2 ■ The distortion-rate function 
prescribes the minimal distortion achievable under the compression rate R12 

min trD. (79) 

D : I{y2;y2)<Ri2 

After evaluating the mutual information expression in ( [79l ) under Gaussian signaling, the minimization in the 
distortion-rate function can be written as 

minimize tr D (80) 
over D G H^' (81) 
subject to logdet(52 - D) < R12. (82) 

Let D* denote the optimal distortion covariance matrix in (ISTl) : it is computed by the reverse waterfilling fTTl 
procedure along the eigenmodes of 52. The 7/2 that achieves the minimum sum distortion in (|80l ) has a joint 
distribution with 2/2 described by 

2/2 = ^2 + 5, y2^CM{0,S2-D*), z^CMiO,D*). (83) 

The joint distribution in ( [83] ) can be equated with the form in (175] ) by setting the corresponding parameters of the 
compression scheme to be 

^RD = (/TV. - D*S^^y/\ Zrd = D* (84) 

where the subscripts in A^d , Zrd are used to designate the compression scheme under consideration. 

C. Wyner-Ziv Compression 

In the compress-and-forward strategy in [4^ Thm. 6], the transmission scheme also takes advantage of the 
correlation between the observed signals at the source and relay using Wyner-Ziv compression. In particular, 
when the destination attempts to reconstruct 1/2 form 1/2, it also has access to its own receive signal yi, which 
can be used to improve the performance of the compression process. The Wyner-Ziv compression approach ll25l 
exploits the correlation between yi and 7/2 as side information at the decoder to achieve a lower compression noise 
level with the same compression rate Ri2. 

With the transmit signal of the source being as specified in (172] ). the covariance of the observed signals at the 
destination (after successive interference cancellation of the relay's signal) and relay, respectively, are given by 

Sn ^ E[ywi\x2] = In, + HuQl^H^^ G H^"^ (85) 

^22 = £[^2^2""] = IN2 + H2iQ*uHi G H^'^ . (86) 
Moreover, the cross-co variance between y2 and yi\x2 is 

^21 = E[y2yf ^2] = H2iQ*uH({ G C^^x^\ (87) 

For Gaussian signals under quadratic distortion, the Wyner-Ziv scheme achieves the same rate-distortion tradeoff 
as if the side information were also present at the encoder 1251 . ||26l (i.e., as if the relay had access to yi in the 
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course of the compression process). Therefore, the Wyner-Ziv compression noise is given by the distortion-rate 
function of compressing the signal 1/2 1 yi using rate R12 

minimize tr D (88) 
over D G H^' (89) 
subject to logdet(S'2|i - D) < Ru (90) 

where D represents the distortion covariance matrix of the Wyner-Ziv compression error, and 521 1 is the conditional 
covariance of 1/2 given yi 

S2\i = ny2y^\yi] = S22 - 8218^,' S^, g h^^^ . (91) 

The parameters for the Wyner-Ziv compression scheme are then identified to be 

Awz = {In. - D*S^^l ) , Z^z = D* (92) 

where D* is the solution to ( [89l ) from reverse waterfilling against S'2|i. Since the side information reduces the 
compression noise, Wyner-Ziv compression always achieves better performance than the rate-distortion compression 
scheme. However, rate-distortion compression has lower implementation complexity, since the correlation between 
the receive signals of the source and relay is not exploited in the compression process. 

D. Numerical Results 

Fig. [To] shows the compress-and-forward rates under the same channel parameters as those considered in Fig. [8] 
The Wyner-Ziv (WZ) compress-and-forward rate outperforms the compression-and-forward rate under rate-distortion 
(RD), which demonstrates the capacity gain from exploiting side information. However, the Wyner-Ziv advantage 
ceases when the relay is close to the destination: in that regime, the efficiency of the compression scheme has 
limited impact, as the relay has a strong channel to the destination. Overall, the compress-and-forward rates do 
not perform as well as the decode-and-forward rates, except when the relay is far from the source and near the 
destination. Moreover, unlike its decode-and-forward counterpart at {dx,dy) = (0,i/io), the compress-and-forward 
rates markedly fall short of the 4x8 MIMO capacity when the relay is at {dx^dy) = (l,Yio). Nevertheless, as 
the relay is under no stipulation to perform any decoding, the compress-and-forward rate is at least as large as that 
under direct transmission, regardless of network geometry. 

VI. Conclusions 

We considered the optimization of transmit signals and bandwidth allocation for MIMO relay channels. We 
assumed that all terminals have channel state information, and we evaluated the cut-set capacity upper bounds 
and the decode-and-forward rates by formulating them as convex optimization problems. The solutions to the 
optimization problems can be efficiently computed by numerical convex optimization methods. In the case of 
half-duplex relaying, where the relay cannot simultaneously transmit and receive in the same frequency band, 
the bandwidth allocation and the transmit signals are jointly optimized. We also presented achievable relaying 
rates based on the compress-and-forward strategy, where the relay does not decode the message from the source, 
but forwards a compressed version of its observation to the destination using the rate-distortion and Wyner-Ziv 
compression schemes. 

When the relay is close to the source, it is observed that the decode-and-forward coding strategy is almost 
optimal: its achieved rate is near the cut-set capacity upper bound, especially in full-duplex relaying. Moreover, 
under the half-duplex constraint, decode-and-forward significantly outperforms orthogonal two-hop relaying. For all 
relaying schemes, the maximum capacity gain over direct transmission is attained when the relay is approximately 
halfway between the source and destination. On the other hand, when the relay is close to the destination, decode- 
and-forward underperforms direct transmission as the source-relay link becomes a bottleneck. In this regime good 
performance is achieved by the compress-and-forward schemes, which always achieve a rate that is equal to or 
better than the direct transmission rate. 
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